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University  of  Texas  at  Austin,  Austin,  Texas 


R.  E.  Olson 

Molecular  Physics  Department 
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Low  energy  (5  -100  eV)  cross  sections  have  been  computed  for  the 

Inelastic  scattering  of  He(2®S)  +  He(l*S)  -  He(2®P)  +  He(l‘s).  Two- 

state  quantum- mechanical  close-coupling  calculations  have  been  performed 

between  the  c®^^  and  b^O  states  of  the  He_  molecule.  The  necessary 
g  g  2 

molecular  wave functions,  potential  curves,  and  rotational  coupling  matrix 
elements  were  obtained  by  ab  initio  computation.  The  inelastic  transi¬ 
tion  probabilities  are  found  to  be  quite  large,  0*®  “ 

impact  parameters  close  to  the  curve-crossing  of  the  c*E^  and  b’ll  poten- 

g  g 

tlal  curves  at  R  »  3.3  a  .  The  inelastic  total  cross-sections  range 
X  o 

from  1.0  X  10"*®cm®  at  5  eV  to  4.0  x  10~*®cm*  at  100  eV.  Reduced  plots 

of  the  inelastic  differential  cross-sections  show  a  sharp  maximum  at  a 

T  =  E0  value  of  approximately  100  eV-deg.  The  strong  inelastic  process 

is  also  found  to  significantly  alter  the  elastic  differential  cross- 

sections  in  the  t  «  25  -125  eV-deg  region.  The  rotational  coupling 

between  the  c®E^  and  b®Tr  states  is  found  to  be  the  dominant  inelastic 
g  g 

process  at  energies  up  to  500  eV. 


^  This  research  was  supported  by  the  Office  of  Naval  Research  under  Con¬ 
tract  No.  N00014-67-012G-0017  and  the  National  Science  Foundation  under 
Contract  GA-33418. 
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Introduction 


Low  energy  (E»5-10  eV)  differential  scattering  cross  sections  of 

nc(2^S)  on  He(l*S)  have  been  measured  experimentally  at  Stanford  Research 

Institute.^  These  cross  sections  show  that  there  exists  an  inelastic  loss 

2 

process  to  the  He(2’P)  +  He(l^S)  channel. 

Analysis  of  the  He,  potentials  shown  In  Fig.  1  indicates  that  there 

are  two  possible  Inelastic  channels  for  loss  of  the  elastic  flux  at  low 

collision  energies.  One  possible  channel  is  to  a  potential  curve  of  the 

state  that  originates  in  He(ls2p  ^P).  This  potential  lies  close  to 
®  rv4- 

the  c^^  potential  curve  over  a  considerable  range  of  Internuclear 

e 

separations.  These  two  states  are  thus  connected  by  radial  (9/dR) 
coupling  matrix  elements.  This  possible  inelastic  loss  process 

3 

has  been  studied  previously  by  Evans,  Cohen  and  Lane.  They  find  that 

the  inelastic  total  cross-sections  Into  the  He(ls2p  *P)  +  He(l^S) 

g 

channel  are  very  small  at  low  enei^les  (Q  =  5.0  X  IC^^cm*  at  5  eV  rising 
to  only  7.7  X  lO"*’  cm*  at  100  eV), 


Another  possible  inelastic  channel  Is  the  b^II  potential  curve 

g  + 

which  also  dissociates  to  He(ls2p  ®P)  +  He(ls‘  ^S)  and  crosses  the  c*2L 

g 

potential  curve  at  approximately  R  =  3.3  a  .  Ibese  two  curves  lie  within 

X  o 

0.5  eV  for  internuclear  separations  of  2.5  a  <R<4.0a.  The  coupling 

o  o 

matrix  element  between  these  states  is  the  rotational  coupling  operator, 
L^.  There  is  a  reasonable  intuitive  basis  for  expecting  this  coupling 
to  be  substantial  since  both  molecular  states  have  3p  molecular  orbital 
configurations.  Thus  this  channel  suggests  itself  as  a  source  for  sub¬ 
stantial  loss  of  flux  from  the  Re(2’S)  -t-  He(l^S)  elastic  channel  at  low 
energies. 


This  paper  presents  close  coupling  computations  on  the  molecular 

channels  c*E^  and  b®n  based  on  ab  liltlo  computed  potential  curves  and 
g  g - 

coupling  matrix  elements.  The  computed  inelastic  total  cross  sections 
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aro  found  to  be  very  large  (0  >  10***  cm*  for  E>5  eV)  and  thus  we 

strongly  suggest  that  this  inelastic  channel  represents  the  major  loss 

2 

process  needed  to  explain  the  experimental  observations. 
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Wave functions.  Potential  Curves  and  Coupling  Matrix  Elements 


The  uave functions  and  potential  cruves  used  to  define  the  molecular 

channels  for  the  cross-section  computations  were  obtained  by  standard 

3  + 

configuration  interaction  methods.  The  function  used  for  the  c  state 
was  a  five  configuration  function  Slater-type  orbital  basis  set  adapted 


from  the  work  of  Greenawalt.  The  wavefunction  is  written  as 


ls*2sl  +  Cjfls^;  ls'2p^J  +  Cjlls^;  ls'3s) 

+  c^[ls2;  ls»3pj  +  Cjtls^;  ls'3d^l, 

where,  in  this  notation  the  orbitals  to  the  left  of  the  semicolon  are  on 

the  ground  state  atom  and  those  to  the  ri^t  the  excited  atom.  The 

2 

brackets  denote  the  formation  of  an  eigenfunction  of  L  ,  S  and  S  .  The 

S  2 

orbital  exponents  were  optimized  at  each  Intemuclear  separation.  The 

potential  curve  obtained  by  these  computations  was  adjusted  by  the 

procedure  of  Klein,  Greenawalt  and  Hatsen.^  The  adjusted  potential  points 

are  given  in  Table  I  along  with  the  potential  points  available  from  the 

more  elaborate  computations  of  Kolker  and  Michels.^ 

3 

The  b  IT  potential  curve  was  computed  with  a  five  configuration 

O 

wavefunction 

fjCbVg)  -  c^[l»2;  ls2p^J  +  Cjtls^;  l83p^l 
+  Cj(ls  ;  Is '3d^)  c^fls  Is';  3d_^J  +  ^2^^* 


Here,  ere  elllpcic  orbitals;  passes  to  Is  in  the  united  atom  limit, 
while  ^2  ^3  3Pq’  orbital  exponents  were  optimized 

at  each  Intemuclear  separation.  The  potential  curve  obtained  by  these 
computations  was  adjusted  by  the  procedure  o£  Klein,  Greenawalt  and  Matsen. 

The  adjusted  potential  curve  is  given  in  Table  I.  The  adjusted  curves  cooipsre 
well  with  the  RKR  curves  of  Ginter  and  Battino.^  A  crossing  point  of 
■>  3.29a^  is  computed  while  Ginter  and  Battino  show  a  crossing  point  of 
approximately  "  3.2  a^.  The  minima  of  both  potential  curves  if  some¬ 
what  larger  than  the  experimental  values.^  The  relative  spacing  seems  to 
be  in  rather  good  agreement.  The  experimental  value  for  the  energy  difference 
between  the  minima  is  0.283315  a.u.  idiile  the  computed  value  is 
0.253315  a.u. 

The  coupling  matrix  elements  <  |  L  I  ^  were  cosiputed 

y  o 

with  Che  standard  package  of  the  Molecular  Physics  Group  at  the  University 
of  Texas  at  Austin.  The  coupling  matrix  elements  were  also  computed  with 
other  rather  different  wavefunctions  for  several  Intemuclear  separations. 

The  values  (listed  in  Table  I)  are  not  sensitive  to  changes  in  the  wave- 
functions. 
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Scattering  Cross  Sections 

In  order  to  calculate  the  cross-sections,  the  potential  points 

given  in  Table  I  were  fit  with  a  cubic  spline  function.  The  coupling 

matrix  elements  were  represented  by  the  exponential  form  =  1.67  x 

exp  (-0.231R}  for  internuclear  separations  R  2  2.77  a^;  within  2.77  a^, 

Ly  was  set  equal  to  0.88.  The  transition  probabilities  and  the  real  and 

imaginary  parts  of  the  S-matrlx  were  then  calculated  at  even  i  (orbital 

angular  momentum  quantum  number)  values  using  the  "amplitude-density" 

8 

method  of  Secrest  and  Johnson. 

The  calculated  transition  probabilities  between  the  and  b’llg 

states  are  shown  in  Fig.  2.  For . convenience  in  comparing  between  ener¬ 
gies,  the  transition  probabilities  are  plotted  versus  impact  parameter, 
b,  where  we  have  used  the  semiclasslcal  relationship  b  =  (jK  -f  l/2)/k,  in 
which  the  wavenumber  k  =  /2hE/A*.  At  5,  8.65,  10,  and  25  eV,  the 
coupled  equations  were  solved  for  the  transition  probabilities  at  all 
even  i  values  so  that  the  real  and  imaginary  parts  of  the  S-matrix 
could  also  be  used  fur  calculating  the  inelastic  differential  cross 
sections.  However,  at  100  eV  the  cost  becomes  prohibitive  so  we  only 
solved  the  coupled  equations  at  Ll  -  50  intervals  from  I  =  100  to 
I  s  700.  These  are  the  points  shown  in  Fig.  2. 

Even  at  the  lowest  collision  energy,  5  eV,  the  transition  proba¬ 
bilities  out  of  the  c^Eg  channel  rise  almost  to  unity.  Furthermore, 
as  the  collision  energy  is  increased,  the  inelastic  process  becomes 
quite  appreciable  at  impact  parameters  larger  than  the  crossing  point 
at  R^  =  3.29aQ.  The  reason  can  be  seen  by  examining  the  dynamical 
nature  of  the  rotational  coupling  matrix  element  that  is  given  in  an 
impact  parameter  formalism  by 

H,,(R)  »  vb  Ly{R)/R> 
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whoro  V  Is  tho  incident  collision  velocity.  From  Eq.  (l)  we  can  see 
that  as  the  energy  Increases  the  coupling  will  become  stronger  and  it 
will  be  possible  to  Induce  Inelastic  transitions  at  larger  and  larger 
Impact  parameters. 

The  transition  probabilities  display  an  unexpected  oscillatory 
dependence  around  b  =  1.5  a^.  There  exists  a  region  of  stationary 
phase  that  causes  a  decrease  in  transition  probability  at  b  1.5  a^  at 
5  oV  and  25  eV  and  added  "shoulders"  at  8.65  and  10  eV.  This  effect  is 
due  to  the  difference  potential  between  the  c’Zg  and  b^IIg  states  pro¬ 
cessing  u  maximum  in  the  region  about  R  =  1.5  a^  (see  Table  l).  We 
would  therefore  expect  to  see  oscillations  on  the  inelastic  total  cross 
sections  due  to  this  region  of  stationary  phase  in  the  phase  differences. 
However,  for  this  case  the  oscillations  will  be  very  weak  because  of  the 
dominance  of  the  large  Impact  parameter  transition  probabilities  in  the 
total  cross  section. 

At  both  25  and  100  eV  there  is  a  large  Increase  in  the  transition 
probabilities  at  small  Impact  parameters,  b  ^  1.0  a^.  This  Increase  is 
due  to  a  second  crossing  of  the  c*£g  and  b^Hg  states  at  around  1.0  a^ 
(Table  l).  The  energy  required  to  reach  this  region  is  approximately 
18  eV  so  the  effect  is  only  observed  at  25  and  100  eV.  Since  the  wave- 
functions  used  for  constructing  the  potential  curves  were  not  set  up  to 
be  extremely  accurate  at  small  internuclear  separations,  there  may  be 
some  slight  changes  in  the  determination  of  the  inner  crossing.  However, 
wc  may  say  with  certainty  that  there  is  an  inner  crossing  which  will 
become  important  at  high  collision  energies,  E  ^  100  eV. 

The  inelastic  tetal  cross  sections  (the  sum  of  both  direct  excita¬ 
tion  and  exchange  excitation)  were  calculated  by  the  summation 

z  (2f  +  l)ls„|®  (2) 

i  even 
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In  Eq.  (2)  we  have  Included  a  statistical  weight  factor  of  one  half 
since  one  half  of  the  Incident  flux  is  going  into  the  nonreactive  un- 
gerade  channel.  At  100  eV  the  Inelastic  total  cross-sections  were 
evaluated  by  interpolating  between  the  calculated  transition  probabili¬ 
ties.  The  calculated  total  cross  sections  are  given  in  Table  II  and  are 

displayed  in  Fig.  3.  Also  shown  in  Fig.  3  are  the  Inelastic  total  cross 

3 

sections  calculated  by  Evans,  Cohen,  and  Lane  for  radial  coupling  to  the 
He(2®P)  +  He(l*S)  state. 

Our  calculations  can  be  extended  to  higher  energies  by  using  the 
10 

method  of  Russek.  Inserting  our  potential  parameters  for  the  outer 
crossing  at  =  3.29aQ,  we  obtain  the  cross  sections  shown  by  the 
dashed  line  in  Fig.  3.  The  close-coupled  values  lie  somewhat  above 
those  evaluated  using  Russek's  model.  However,  this  is  expected  since 
the  crossing  potentials  are  not  linear  in  the  region  around  R^,  as 
assumed  in  the  parametrization  of  the  model,  but  actually  lie  very  close 
to  one  another,  Fig.  1  and  Table  I.  We  therefore  expect  Increased 
coupling  and  hence  larger  transition  probabilities  and  inelastic  total 
cross-sections  for  our  case. 

There  also  appears  to  be  a  hint  of  oscillatory  structure  on  our 
cross  sections,  which  correlates  to  the  stationary  phase  region  of  the 
transition  probabilities.  The  8.65  and  10  eV  cross  sections  lie  some¬ 
what  farther  above  the  dashed  line  than  the  5  and  25  eV  calculations 
indicating  a  possible  maximum  in  the  energy  dependence  of  the  cross 
sections.  The  100  eV  cross  section  deviates  significantly  above  the 
estimates  using  Russek's  method  and  this  is  due  to  added  contribution 
to  the  transition  probabilities  from  the  inner  crossing. 

If  we  use  Russek's  method  as  a  guide  for  extrapolating  our  inelas¬ 
tic  total  cross-sections  to  higher  energies,  we  obtain  a  value  of 
Q  ^  6  X  10**^°  cm^  at  1000  eV.  Comparing  our  results  to  the  calculations 
of  Ref.  3,  where  radial  coupling  was  considered  as  Inducing  the 
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inolasllc  transitions,  it  is  evident  that  rotational  coupling  will  be 
the  dominant  inelastic  process  up  to  collision  energies  of  at  least 
500  eV. 


The  inelastic  direct  excitation  differential  cross  sections  were 
evaluated  from  the  real  and  imaginary  parts  of  the  S-matrlx  using  the 
formulas 

0^3(0)  =  (A*  +B»)/(4k»)  (3) 

where 


A  =  S 
even 


(2Jfc  +  1)  P/cose) 


(4) 


and 

B  =  I  (2f  +  1)  ImCs^g)  P^(cose)  .  (5) 

even  i 

The  results  of  these  calculations  are  shown  in  Fig.  4.  The  inelastic 
differential  cross-sections  are  plotted  to  90^  and  not  to  180**  since 
they  are  symmetric  about  90°.  The  major  characteristics  of  these  cross- 
sections  are  a  peak  of  large  magnitude  at  small  angles  with  super¬ 
imposed  rainbow  structure  and  the  prescence  of  nuclear  symmetry  oscilla¬ 
tions  at  large  angles. 

A  more  eonvenient  way  to  plot  the  inelastic  differential  cross- 
section  is  to  use  the  reduced  coordinates  t  =  E0  and  p(0)  =  9  sin  0  o(0). 
The  rapid  decrease  in  the  cross-section  is  removed  by  using  the  reduced 
cross-section  p.  The  reduced  angle  t  Is  convenient  because  it  is 
inversely  related  to  the  Impact  parameter  (large  T,  low  b;  etc.).  Utili¬ 
zing  this  type  of  plot,  there  are  two  major  peaks  found  in  the  forward 
scattering  cross  sections  (O®  <  0  <  90®).  The  largest  peak  is  centered 
at  about  100  eV-degrees,  Fig.  5,  and  corresponds  to  the  transition 
probabilities  found  for  impact  parameters  in  the  range  of  about  b  pa  3.0  a^. 
Fig.  2.  The  second  peak  (centered  at  T  *  175  eV-deg  at  5  eV,  160  eV-deg 
at  8.65  cV,  and  150  eV-deg  at  10  eV)  is  a  Stuckelberg  oscillation  due 
to  an  interference  caused  by  the  two  possible  trajectories  within  R^. 
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Superimposed  on  the  cross  sections  displayed  in  Fig.  5  are  fine 
structure  oscillations  due  to  inelastic  rainbow  scattering  out  of  the 
deep  potential  wellj  e  2.2  eV  (Table  l).  This  structure  is 

prevalent  only  at  small  angles,  300  eV-deg,  and  at  the  larger  angles 
nuclear  symmetry  oscillations  are  observed.  Fig.  4.  The  nuclear  symme¬ 
try  oscillations  are  caused  by  the  Indistinguishabillty  of  the  colliding 
particles  and  the  fact  that  it  is  impossible  to  differentiate  between 
direct  excitation  scattering  at  angle  6  and  excitation  transfer  from 
n-6.  The  nuclear  symmetry  oscillations  have  their  greatest  magnitude 
at  90°.  Similar  sets  of  oscillations  with  the  same  interpretation  also 

occur  on  low  energy  He"*^  +  He(l*s)  -♦  He'*’  +  He(2®S)  inelastic  differential 
11 

cross  sections. 

Another  interesting  phenomenon  displayed  on  the  cross  sections  of 

Fig.  5  is  the  energy  dependence  of  the  magnitude  of  the  inelastic 

rainbow  oscillations.  The  magnitude  of  the  fine  structure  on  the  5  eV 

and  especially  on  the  25  eV  cross  sections  is  much  less  than  that  of 

the  8.65  and  10  eV  cross  sections.  These  magnitude  fluctuations  are 

caused  by  the  region  of  stationary  phase  at  b  1.5  a^  in  the  transition 

probabilities.  Fig.  2.  The  inelastic  rainbow  oscillations  require 

inelastic  scattering  on  the  positive  angle  branch  of  the  special  deflec- 
12 

tion  functions  for  interference  effects  to  occur  between  the  trajec¬ 
tories  scattering  to  a  negative  angle  and  a  positive  value  of  the  same 
angle.  For  this  system  the  positive  angle  scattering  arises  from 
Impact  parameters  b  1.7  a^.  Therefore,  the  magnitude  of  the  transi¬ 
tion  probabilities  in  the  b  ss  1.7  a^  region,  Fig.  2,  can  be  directly 
correlated  to  the  magnitude  of  the  fine  structure  oscillations. 

Elastic  2-state  and  3-state  differential  cross  section  calculations 
wore  also  performed  to  ascertain  the  effect  of  the  inelastic  channel  on 
the  elastic  cross  sections.  In  this  calculation  we  employed  an  a^Ej^ 

potential  that  was  compatible  with  rainbow  scattering  measurements^, 

13  4. 

the  latest  spectroscopic  analysis  on  the  A*Zu  state,  and  the 

10 


J.*!  4 

calculations  of  Liu  on  the  well  depth  oi'  Ilo.j  ,  The  cross  sociioii 
formulas  have  been  given  previously. 

In  l-'lg.  6  arc  shown  the  results  of  the  clastic  calculations  with 
(3-slalt»)  and  without  (2-stato)  the  inclusion  of  the  inelastic 

channel.  There  is  very  little  change  in  the  position  of  the  main 
rainbow  from  the  state  at  t  w  160  eV-deg,  At  smaller  angles,  how¬ 

ever,  25  ^  T  125  oV-dog,  the  inelastic  channel  considerably  modifies 
the  cross  sections.  This  is  due  to  the  loss  of  elastic  flux  out  of  the 
b^Eg  state  and  the  "perturbation"  structure  caused  by  the  inelastic 
transition. 
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Conclusions 


In  conclusion,  we  have  calculated  the  ^  init  io  pott»iu  tal 
curves  for  the  c’Zg  state  of  He(2®S)  +  He(l's)  and  the  b’Hg  state 
of  llo(2®P)  +  IIc(l*s)  alonp:  with  the  rotational  coupling  matrix 
clement,  L  .  It  is  found  that  the  inelastic  transition  probabilities 

y 

induced  by  rotational  coupling  between  the  c®Eg  and  b^IIg  molecular 

states  arc  quite  large  even  at  low  energies,  5  eV,  and  dominate  the 

inelastic  scattering  up  to  at  least  500  eV. 

The  inelastic  total  cross  sections  (the  Siam  of  both  direct  and 

exchange  excitation)  were  calculated  and  found  to  range  from 

1.0  X  10“*®  cm^  at  5  eV  to  4.0  x  10  *®  cm*  at  100  eV.  The  inelastic 

differential  cross-sections  were  found  to  peak  at  a  reduced  angle  t 

of  about  100  eV-deg.  This  has  been  conflirmed  by  a  recent  experiment 

2 

performed  at  8.65  eV.  Nuclear  symmetry  oscillations  were  superimposed 
on  the  inelastic  differential  cross-sections  at  large  angles*  At  small 
angles,  T  ~  300  eV-deg,  inelastic  rainbow  structure  was  observed.  The 
clastic  differential  cross  sections  were  also  calculated  and  were  found 
to  be  significantly  altered  by  the  inelastic  channel  in  the  angular 
region  25  ~  t  .s,  125  eV-deg, 
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Table  1 


Potential  Energy  Curves  and  Coupling  Matrix  Elements 


<  1  L  1  ■ 

g  y 

1.0 

-4.40otr43 

-4.408026 

0.937885 

1.5 

-4.989949 

-5.027520 

0.898538 

1.75 

-5.063993 

-5.090993 

0.893320 

2.0 

-5.089264 

-5.115047 

0.890965 

2.5 

-5.085316 

-5.103260 

0.879553 

3.0 

-5.070398 

-5.05483 

-5.078246 

0.839906 

3.25 

-5.067055 

- 

-5.067487 

0.799513 

3.5 

-5.065814 

-5.05483 

-5.058907 

0.750614 

3.75 

-5.066228 

-5.052358 

0.706118 

4.0 

-5.067162 

-5.05943 

-5.047515 

0.664079 

5.0 

-5.072013 

-5.06553 

-5.038799 

0.549118 

6.0 

-5.075204 

-5.06913 

-5.036836 

8.0 

-5.036368 

10.0 

-5.078784 

-5.07223 

-5.036411 

15.0 

-5.078949 

-5.035280 

a)  This  calculation  adjusted  by  the  procedure  of  ref.  5 

b)  text  ref.  6. 
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Table  II 

INELASTIC  TOTAL  CROSS  SECTIONS 


Energy -eV 


Q  -  10~**  cm* 


5 

1.04 

8.65 

1.42 

10 

1.45 

25 

2.05 

100 

4.0 
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FIGURE  CAPTI(»fS 


Potential  curves  from  Glnter  and  Battlno  (Ref.  7/  leading  to 

the  separated  atom  limits  of  He(2’s)  +  He(l*S)  and  He(2’P)  +  He(l*S). 

Transition  probabilities  vs.  Impact  parameter  for  the  collision 
energies  of  5,  8.65,  10,  25,  and  100  eV.  At  100  eV  the  calculations 
were  performed  at  a  wide  grid  of  impact  parameters,  these  points 
arc  depicted  by  the  solid  circles.  To  calculate  the  total  inelastic 
cross  section,  the  100  eV  transition  probabilities  were  crudely 
interpolated.  This  uncertainty  is  Illustrated  by  the  dashed  line. 

Calculated  rotational  Inelastic  total  cross  sections;  solid  circles. 
The  calculations  of  Evans,  Cohen,  and  Lane  (Ref.  3)  for  radial 
coupling  arc  shown  by  the  solid  squares.  An  estimate  of  the  rota¬ 
tional  Inelastic  total  cross  sections  using  the  method  of  Russ^ 

(Ref.  10)  and  our  potential  parameters  (Table  I)  is  given  by  the 
dashed  line. 

Inelastic  differential  cross  sections  as  a  function  of  angle  for 
the  collision  energies  of  5,  8.65,  10,  and  25  eV. 

Reduced  inelastic  differential  cross  sections  p  «  9  slnQc78  as 
a  function  of  reduced  angle  t  =  Ed  for  the  collision  energies  of 
5,  8.65,  10,  and  25  eV. 

Reduced  elastic  differential  cross  sections  p  as  a  function  of  reduced 
angle  T.  The  dashed  line  is  a  2-state  calculation  using  the  adiabatic 
potentials  a^E'*’  and  that  lead  to  He(2’S)  +  He(I*S).  The  solid 

line  is  a  S-stHte  calculation  that  includes  the  b’ll  inelastic 
channel.  ® 


Figure  2 


0.20 


Figure  5, 


